T  pOCUMENTATION  PAGE 


2b.  DECLASSIFICATION  /  OOWNG 


4.  PERFORMING  ORGANIZATION  REPORT  NUM8£R(S 


6a.  NAME  OF  PERFORMING  ORGANIZATION  1 6b.  OFFICE  SYMBOL 

I  (If  tppikible) 

Iowa  State  University  I 


6c  ADORESS  (Cty,  Sure,  end  ZIP  Code) 

Dept,  of  Aerospace  Engineering 
Ames,  IA  50011 


3*.  NAME  OF  FUNOING/  SPONSORING 
ORGANIZATION 

AFOSR 


8c  ADORESS  (City,  SUt t.  tod  ZIP  Code) 


8b.  OFFICE  SYMBOL 
(If  tppiksbh) 

NA 


BUILDING  410 

BOLLING  AFB ,  DC  20332-6449 


1 1  TITLE  (Include  Security  Geatfketion) 


1b.  RESTRICTIVE  MARKINGS 

-  V  l  ■»,  •  *' 


3  OiSTRiaUfiON7AVAiLA8IUTYOfSlPoST^^^^^^ 

APPROVED  FOR  PUBLIC  RELEASE 
DISTRIBUTION  IS  UNLIMITED 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

_ *FOSK • {*■  n  9  -  u  a  a 7 


7a.  NAME  OF  MONITORING  ORGANIZATION 

AFOSR/NA 


7b  ADORESS  (Cty,  Stste.  tod  ZIP  Code) 

BUILDING  410 

BOLLING  AFB f  DC  20332-6448 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

AFOSR  Grant  85-0357 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM  PROJECT 

ELEMENT  NO.  NO. 

61102F  2-207 


"Three  Dimensional  High  Speed  Boundary  Layer  Flows" 


12.  PERSONAL  AUTHOR(S) 

George  R.  Inger 


13a.  TYPE  OF  REPORT 


16.  SUPPLEMENTARY  NOTATION 


C0SAT1  COOES 


FIELD  I  GROUP  I  SUB-GROUP 


18  SUBJECT  TERMS  (Continue  on  reverse  if  necesxery  end  identify  by  block  number) 

Viscous-Inviscid  Interactions 
High  Speed  Turbulent  Boundary  Layers 
Vortex  Arrays  >  .  - — 


19.  ABSTRACT  (Continue  on  reverse  if  necesury  end  identify  by  block  number) 


This  final  report  summarizes  the  research  results  obtained  by  a  two-faceted  basic 
theoretical  investigation  involving  11)  fundamental  analyses  of  three-dimensional 
viscous- inviscid  interaction  effects  within  high  speed  turbulent  boundary  layers  and 
(2)  study  of  the  influence  of  streanwi.se  vortex  arrays  on  the  skin  friction  within 
attached  or  separated  laminar  layer  fLows.  The  work  described  was  carried  out  from 
September  1985  to  May  31,  1988. 


20.  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT  21  ABSTRACT  SECURITY  CLASSIFICATION 

□  UNCUSSIFIEDAJNLIMITED  SAME  AS  RPT  □  OTIC  USERS  UNCLASSIFIED 


22*.  NAME  OF  RESPONSIBLE  INOIVIOUAL  22b  TELEPHONE  (Include  tree  Code)  22c.  OFFICE  SYMBOL  I 

LEN  SAKELL  202-767-4935  AFOSR/NA] 


00  FORM  1471, 14  MAR  83  APR  edition  may  be  uied  until  exhausted. 

All  other  edition*  are  obsolete. 


H  avuuar  *->  > ;  i  ix* 


TION  OF  ’HIS  PAG 


A>  Oi>K  •  7K . 


&  9  ~  u  28? 


AFOSR  GRANT  85-0357 

"Three  Dimensional  High  Speed  Boundary  Layer  Flows" 

FINAL  TECHNICAL  REPORT 

G.  R.  Inger 
Iowa  State  University 
Ames,  Iowa  50011 

September  28,  1988 


TABLE  OF  CONTENTS 


Page 


ABSTRACT 

INTRODUCTION .  1 

SUMMARY  OF  ACCOMPLISHMENTS 

2.1  Studies  of  Boundary  Layer  Vortex  Effects .  1 

2.2  Theoretical  Investigation  of  Three-Dimensional 

Interactions .  2 

2.3  Complementary  Activi ties .  6 

SUMMARY  AND  CONCLUSIONS .  8 

REFERENCES .  9 

FIGURES . 11 

APPENDICES 


Accession  For 

_  / 

NTIS  GEA&I 

1ST 

DTIC  IA3 

r—i 

Uu  ‘iiu ; :  >  ir;  o.  e  d 

n  1 

J  '0 ^  \  I  0  '  1  . _ 

— j 

A'1 


1. 


ABSTRACT 

This  final  report  summarizes  the  research  results  obtained  by  a  two-faceted 
basic  theoretical  investigation  involving  (1)  fundamental  analyses  of  three- 
dimensional  viscous-inviscid  interaction  effects  within  high  speed  turbulent 
boundary  layers  and  (2)  study  of  the  influence  of  streamwise  vortex  arrays 
on  the  skin  friction  within  attached  or  separated  laminar  layer  flows.  The 
work  described  was  carried  out  from  September  1985  to  May  31,  1988. 
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1.  INTRODUCTION 

The  overall  objective  of  this  research  was  the  basic  theoretical  investigation 
of  three-dimensional  pressure,  skin  friction,  and  flow  field  disturbances 
in  both  laminar  and  turbulent  boundary  layer  flows  including  viscous-inviscid 
interaction  effects,  separation  and  reattachment.  A  sound  understanding 
of  these  phenomena  is  required  in  modern  aerodynamic  design  analyses  of  high-speed 
flight  vehicles.  Indeed,  the  interest  in  and  need  for  such  research  appears 
to  be  even  stronger  now  than  when  the  current  investigation  was  first  begun. 

Broadly  speaking,  the  primary  emphasis  in  our  studies  has  been  to  seek 
a  basic  physical  understanding  of  the  underlying  fluid  behavior  by  means 
of  analytically-oriented  methods;  in  this  way,  the  results  can  be  used  to 
guide  and  interpret  concurrent  experimental  and  computationally  oriented 
investigations.  Specifically,  our  inquiry  has  focused  on  two  parallel  paths 
of  investigations;  (1)  three-dimensional  viscous-inviscid  interaction  phenomena 
within  turbulent  boundary  layers  in  supersonic  flow  due  to  impinging  swept 
shock  and/or  3-D  surface  deflections;  and  (2)  streamwise  vortex-disturbance 
mechanisms  within  boundary  layers  that  are  either  separated  or  attached. 

The  project  has  proven  very  successful  both  in  solving  a  number  of  the  target 
flow  problems  and  in  developing  new  ideas  and  tools;  this  report  summarizes 
these  results. 
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2.  SUMMARY  OF  ACCOMPLISHMENTS 
2.1.  Studies  of  Boundary  Layer  Vortex  Effects 

A  major  phase  of  the  streamwise  vortex  investigation  culminated  in  applying 

1-5 

our  previously-developed  theory  of  disturbance  vortex-pair  formation  in 
separating  boundary  layer  flows  to  the  analysis  of  streamwise  vortex  structure 
in  the  near  wake  of  stalling  wings  operating  at  low  Reynolds  numbers  (see 
Fig.  1).  A  paper  on  this  work  was  presented  to  the  Royal  Aeronautical  Society 
at  its  October  1986  meeting  on  Low  Reynolds  Number  Aerodynamics  in  London 
and  has  been  published  in  the  Proceedings  of  the  Meeting.** 

A  second  achievement  was  the  completion  of  a  detailed  analytical  study 
of  spanwise,  periodic  3-D  disturbances  (representative  of  an  array  of  counter-rotating 
pairs  of  streamwise  vortices)  within  an  entire  family  of  host  pressure-gradient 
laminar  boundary  layer  flows,  namely  the  self-similar  incompressible  Faul kner-Skan 
flows  along  a  wedge  (Figure  2).  In  addition  to  correcting  and  explaining 
several  key  features  of  Fannelop's  earlier  study  of  the  flat  plate  limit 
of  this  problem, ^  our  new  results  significantly  extend  the  basic  understanding 
of  how  streamwise  vortices  alter  the  flow,  skin  friction,  and  incipient  separation 
behavior  of  nominally  2-D  boundary  layer  flows.  This  work  was  documented 
by  the  M.S.  thesis  of  Mr.  M.  Konno  and  has  been  submitted  for  publication 
(an  abstract  of  this  paper  is  included  as  Appendix  A  of  this  report). 

The  aforementioned  work  was  subsequently  extended  to  treat  the  situation 
where  the  streamwise  vortices  are  generated  at  the  wall  by  small  spanwise 
ripples  in  its  surface  that  grow  linearly  in  amplitude  in  the  downstram  direction. 

The  results  were  presented  in  a  paper  at  the  Summer  1986  U.S.  National  Fluid 
Mechanics  Conference,  a  copy  of  which  is  attached  as  Appendix  B. 
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2.2  Theoretical  Investigation  of  Three-Dimensional  Turbulent  Interact ions 


This  aspect  of  our  research  had  three  major  goals: 

•  Apply  modern  asymptotic  methods  (high-Re)  to  elucidate  the  basic  flow 
field  physics  of  3-D  interactions,  especially  near  the  surface; 

•  Extract  the  fundamental  similtude  laws  governing  swept  interactions 
within  the  framework  of  "Law  of  the  Wall /Wake"  concepts  that  characterize 
the  incoming  turbulent  boundary  layer; 

•  Assess  the  validity  of  the  "Independence  Principle"  approximation 
for  some  representative  compressible  turbulent  interactive  flows. 

In  general  support  of  these  studies  we  first  completed  a  basic  theoretical 

study  of  nonisotropic  turbulence  model  effects  on  the  Law  of  the  Wall  region 

of  nonseparating  three-dimensional  turbulent  boundary  layers  (this  work  was 

subsequently  factored  into  our  analyses  of  swept  shock/turbulent  boundary 
8  9 

layer  interaction  ’  ).  A  paper  on  this  work  was  prepared  for  the  January 
1987  AIAA  Aerospace  Sciences  Meeting  (Fig.  3). 

Next,  a  major  emphasis  was  placed  in  recent  years  on  further  refinement 

O 

of  our  original  swept-shock  interaction  theory  for  cyl indri call y- symmetric 
flows  (Fig.  4).  Our  approach  is  based  on  the  3-D  extension  of  our  previously 
successful  nonasymptotic  triple-deck  interaction  theory  (Fig.  5),  combined 
with  the  aforementioned  nonisotropic  3-D  viscosity  model.  A  paper  on  the 
application  of  this  theory  to  the  swept  compression  ramp  problem  was  presented 
in  July  1985  to  the  AIAA  18th  Fluid  and  Plasma  Dynamics  Conference^  and 
subsequently  in  September  1985  to  the  IUTAM  Symposium  on  Turbulent  Shear 
Layer/Shock  Wave  Interactions  held  in  Paris,  being  published  in  the  proceedings 
of  the  latter.^  Among  the  many  new  analytical  results  obtained  were  (a) 
relationships  governing  the  generation  of  significant  secondary  flow  vorticity 
by  the  lateral  pressure  gradient  and  boundary  layer  cross  flow  effects 
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(Fig.  6);  (b)  assessment  of  the  validity  of  the  independence  principle  and 

the  influence  of  deviations  from  it  due  to  compressibility  and  cross-flow 

turbulence  effects  as  a  function  of  sweep  angle,  yielding  predictions  in 

good  agreement  with  experiment  (Fig.  7);  (c)  examination  of  the  3-D  skin 

friction  line  behavior  and  approach  to  separation  (Fig.  8),  including  a  satisfactory 

prediction  of  the  enhanced  incipient  separation  due  to  shock  sweep  (Fig. 

9);  and  (d)  a  new  universal  solution  for  the  inner  viscous/turbulent  shear 
disturbance  deck  that  accounts  exactly  for  the  cross-flow  to  streamwise  eddy 
viscosity  ratio  effect  over  the  entire  range  of  sweep  angles.  The  similtude 
properties  of  these  swept  interactions  as  regards  Mach  Number  and  Reynolds 
Number  Effects  were  also  investigated,  reported  on  at  a  1986  AIAA  Meeting 
and  documented  in  AIAA  Paper  86-0395  (Ref.  11). 

Our  investigation  was  next  broadened  to  consider  the  case  of  coni  cal ly-symmetri c 
shock  interactions  (Fig.  10).  Since  it  is  an  important  experimental  question 
and  of  basic  theoretical  significance  as  well,  we  particularly  addressed 
the  problem  of  the  so-called  "inception  distance"  in  such  interactions  (see 
Figure  11).  A  fundamental  analysis  of  the  governing  vi scous-invi scid  interaction 
equations  in  an  appropriate  shock-oriented  spherical  coordinate  system  revealed 
a  number  of  important  physical  features  concerning  the  radial  and  far-field 
behavior  of  such  interactions.  Moreover,  this  analysis  served  to  derive 
an  experimentally  verified  prediction  of  the  inception-length  dependence 
on  the  upstream  influence  scale  normal  to  the  shock  and  the  sweep  angle. 

A  paper  on  this  work  was  presented  at  the  1986  USAF/FRG  D.E.A.  meeting  at 
Dayton,  Ohio  and  the  10th  U.S.  National  Applied  Mechanics  Congress  at  the 
University  of  Texas  (both  in  June  1986)  and  subsequently  published  in  early 
1987  by  the  AIAA  Journal.12 
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The  aforementioned  study  of  the  outboard  radial -inception  length  problem 
for  swept  shocks  has  been  accompanied  by  an  analytical  study  focused  on  the 
inboard,  highly-nonconical  (or  noncyl indrical ly-symmetric)  interactive  region 
near  the  shock  generator  in  order  to  gain  a  basic  understanding  of  its  interactive 
physics  and  the  resulting  inception  distance  phenomena:  as  shown  clearly 
by  our  far  field  analysis,  only  detailed  study  of  this  region  can  reveal 
the  specific  quantative  dependence  of  ^ncept  on  Mach  Number,  Reynolds  Number, 
and  shape  of  the  incoming  boundary  layer  profile.  To  fix  ideas,  this  study 
addressed  the  specific  problem  of  the  corner  region  formed  by  a  thin  fin 
on  a  wall  (Figure  12)  for  a  supersonic,  adiabatic  laminar  unseparated  interaction. 
The  study  was  being  carried  out  as  the  Ph.D.  dissertation  work  of  Mr.  S.  Ahmed 
using  a  3-D  triple-deck  and  double  Fourier-transform  method  of  approach. 

At  the  time  of  contractual  termination,  the  entire  spectrum  of  the  solution 
in  the  Fourier  Plane  had  been  completed  and  was  undergoing  inversion  to  obtain 
the  streamwise  and  lateral  physics  including  the  upstream  influence  characteristics 
and  skin  friction  lines. 

Further  understanding  of  the  underlying  physical  trends  governing  swept 

shock-boundary  layer  interactions  was  obtained  by  working  out  a  direct  comparison 

of  our  turbulent  triple  deck  theory  with  the  M*,  -  2.85  small  compression 

13 

corner  interaction  2-D  flow  data  of  Settles.  As  shown  in  Figure  13,  the 
theory  is  in  excellent  agreement  with  experiment  over  a  wide  range  of  large 
Reynolds  Numbers.  Also  in  evidence  in  this  Figure  is  the  fact  that  such 
comparisons  are  sensitive  to  the  shape  of  the  incoming  turbulent  boundary 
layer  profile.  Accordingly,  we  have  also  made  a  detailed  study  of  thi s  question 
by  recasting  our  entire  theory  in  terms  of  Law  of  the  Wall/Law  of  the  Wake 
Parameters  and  then  examining  the  influence  of  the  Outer  Wake  Function  and 
Inner  Deck  Turbulence  Modeling,  respectively,  on  the  predicted  interaction 
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properties.  Typical  results  taken  from  the  documenting  AIAA  Paper  (attached 
as  Appendix  C),  are  illustrated  here  in  Figures  14  and  15  for  the  upstream 
influence  distance  of  2-D  interactions;  these  results  clearly  indicate  that 
the  wake  function  in  particular  must  be  accurately  known  to  affect  meaningful 
experimental/theoretical  comparisons  in  high  Reynolds  Number  flows.  This 
work  has  had  a  significant  impact  on  the  planning  of  future  experimental 
and  CFD  studies  of  shock/boundary  layer  interactions  (see  Appendix  D). 

The  aforementioned  Law  of  the  Wall/Wake  similtude  version  of  our  turbulent 
triple  deck  theory  was  subsequently  extended  to  3-D  swept  interactions  with 
either  cylindrical  or  conical  symmetry.  Typical  results  from  this  work, 
which  was  documented  in  Ref.  14,  are  given  in  Figures  16  and  17  to  illustrate 
the  influence  of  Mach  number  and  shape  factor,  respectively,  over  the  sweepback 
effect  on  interactive  upstream  influence  distance. 

2.3  Complementary  Activities 

As  an  outgrowth  of  an  earlier  AFOSR-sponsored  two-month  visit  in  1984 
by  Dr.  Nandanan  (a  previous  collaborator),  in  which  we  devised  a  new  theoretical 
approach  for  predicting  interaction-zone  slot  suction  effects  on  supercritical 
airfoils,  a  joint  paper  was  prepared  on  our  work  and  subsequently  presented 
to  the  May  1986  AIAA/ASME  Fluid  Mechanics  and  Plasma  Dynamics  Meeting  in 
Atlanta  (Figure  18).  Dr.  Nandanan  was  also  invited  to  present  the  work  at 
the  Summer  of  1986  Asian  Fluid  Mechanics  Meeting  in  Tokyo.  This  exploratory 
research  could  prove  a  significant  guide  to  future  studies  of  the  use  of 
localized  suction  control  of  3-D  interactions. 

Another  auxiliary  investigation  completed  under  the  partial  auspices 
of  this  contract  was  an  extension  of  the  triple  deck  theory  to  treat  moderately 
hypersonic  shock/laminar  boundary  layer  interactions.  The  results,  which 
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were  presented  to  an  AIAA  meeting  and  documented  in  Paper  88-0603  (Figure  19), 
will  serve  to  guide  subsequent  studies  of  hypersonic  turbulent  3-D  interactions. 

Finally,  it  is  important  to  mention  a  significant  and  unique  aspect 
of  this  research  throughout  its  execution:  the  AFOSR-encouraged  cooperation 
with  selected  other  investigators  working  in  parallel  on  the  problem,  namely 
Professors  G.  Settles,  D.  Dolling  and  S.  Bogdonoff  (experimental  studies) 
and  Dr.  C.  M.  Horstman  and  Prof.  D.  Knight  (CFD  studies).  This  further  resulted 
in  a  newly-formed  AIAA  Ad-hoc  Hirh  Speed  Flow  Panel  that  meets  twice  a  year 
at  the  Summer  and  Winter  AIAA  Meetings,  and  also  in  two  successive  two-day 
summer  3-D  Shock/Turbulent  Boundary  Layer  Interaction  Workshops  at  Penn  State 
(1987)  and  Princeton  University  (1988). 
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3.SUH<ARY  AND  CONCLUSION 


The  triple-deck  approach  embodied  in  this  interaction  study  part  of 
this  project  has  enabled  the  successful  analysis  of  non-separating  3-D  SBLI's 
including  the  basic  features  of: 

•  Turbulence  modeling  in  the  middle  deck  (wake  function  effect); 

•  Turbulence  modeling  in  the  inner  deck  (Law  of  the  Wall  aspects,  including 
non-isotropy  effects  and  the  influence  of  the  cross-flow  turbulence); 

t  Streamwise  vorticity  generation  and  skin  friction  line  geometry 
prediction; 

t  Determination  of  the  similtude  laws  for  §Q  &  Cp^-  s  = 

FnS(a»  and  tt ^ )  over  a  wide  range  of  conditions. 

In  particular,  our  results  have  yielded  fundamental  understanding  of  sweep 
effects,  showing  that:  (a)  upstream  influence  increases  rapidly  with  sweepback, 
a  major  contributor  to  this  being  the  cross-flow  interactive  turbulence  effect; 
(b)  significant  deviations  from  independence  principle-predictions  can  occur 
for  A  >  30°,  due  to  a  combination  of  compressibility  effects  on  .ariable 
fluid  properties  and  the  interactive  cross-flow  turbulence  effect;  (c)  a 
pronounced  decrease  in  incipient  separation  pressure  (i.e.  "earlier"  separation) 
with  sweep  is  predicted  in  good  agreement  with  experiment. 
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ABSTRACT 

"THEORY  OF  STEADY  SPANWISE-PERIODIC 
DISTURBANCES  IN  FALKNER-SKAN  BOUNDARY  LAYERS" 

G.  R.  Inger  and  M.  Konno* 

Iowa  State  University 
Ames,  Iowa  50011 
United  States  of  America 

This  work  describes  a  theoretical  investigation  of  the  effect  of  weak 
spanwise-per iodic  cross  flow  disturbances  in  the  freestream  above  a  high 
Reynolds  number  laminar  boundary  layer  on  a  wedge-shaped  body.  Our  emphasis 
is  on  a  careful  treatment  of  the  spanwise-gradient  viscous  effects  in  the 
perturbation  field  for  arbitrary  values  of  the  spanwise  wavelength.  In 
addition  to  explaining  and  eliminating  a  heretofore-puzzling  singularity 
encountered  in  previous  studies  of  this  flow  problem,  we  also  bring  out 
how  the  host  flow  pressure  gradient  influences  the  3-D  disturbance  effect 
on  such  properties  as  the  skin  friction  and  displacement  thickness  distributio 


: k 

Glenn  Murphy  Distinguished  Professor  and  Graduate  Research  Assistant, 
respectively,  Dept,  of  Aerospace  Engineering. 


EXTENDED  SUMMARY 


"THEORY  OF  STEADY  SPANWISE- PERIODIC 
DISTURBANCES  IN  FALKNER-SKAN  BOUNDARY  LAYERS" 

1 .  Introduce  ion 

Thorough  understanding  of  3-D  disturbance  effects  on  a  laminar  boundary- 
layer  is  an  important  problem  in  fluid  mechanics.  In  particular,  it  is 
of  interest  in  the  aerodynamic  design  of  both  aircraft  and  turboraachinery 
to  formulate  an  analytical  model  of  streamwise  vorticity  effects  on  the 
properties  of  an  underlying  boundary  layer. 

Crow  [1]  investigated  a  weak  periodic  transverse  flow  effect  on  a 
flat-plate  boundary  layer  using  a  perturbation  expansion  of  the  full  Navier-Stok 
equations  for  steady,  incompressible,  viscous  flow.  Fannelop  [2]  then 
simplified  the  governing  equations  for  this  problem  by  a  boundary- layer 
type  approximation  (i.e.,  O  and^e^  W/tf’x  *-  —  0)  and  obtained  an  equivalent 

result  as  Crow.  In  both  works,  the  spanwise  viscous  derivatives  and  the 
perturbation  pressure  were  neglected,  and  the  resulting  predictions  of 
three-dimensional  effects  were  shown  to  grow  unbounded  downstream,  suggesting 
a  breakdown  of  their  analysis. 

The  present  work  analyzes  the  effects  of  a  weak  steady  periodic  transverse 
flow  on  the  incompressible  laminar  Falkner-Skan  boundary  layer  along  a 
wedge  [Fig.  1].  We  formulate  the  problem  by  also  adopting  a  high  Reynolds 
number  boundary  layer  approach  and  introducing  a  small  perturbation  expansion 
in  the  simplified  governing  equations.  However,  in  our  approach  we  retain 


Figure  1:  Sketch  of  the  flow  model  and 
coordinate  system. 


and  carefully  treat  the  pressure  perturbation  and  the  viscous  z-derivative 
effects  for  arbitrary  spanwise  wave  lengths,  since  the  downstream  disturbance 
development  strongly  depends  on  these  effects  as  shown  later. 

We  first  re- inves t igate  the  limiting  case  of  flatplate  flow  to  explain 
and  eliminate  the  solution  singularity.  Correct  numerical  solutions  are 
obtained  using  an  iterative  outward  shooting  method  developed  by  Nachtsheim 
and  Swigert  [3].  Qualitative  behavior  of  the  resulting  key  boundary  layer 
properties  including  the  skin  friction  and  three-dimensional  displacement 
boundary  layer  thickness  are  then  presented.  The  work  is  then  further 
generalized  to  study  the  host  flow  pressure  gradient  effect  on  the  disturbed 
boundary  layer  properties.  Well-behaved  solutions  are  obtained,  and  the 
disturbance  velocities  and  shear  stresses  under  both  favorable  or  adverse 
pressure  gradient  of  different  intensities  are  presented  and  discussed. 

2.  Formulation  of  the  Analvsis 


The  following  analysis  considers  a  nominally  two-dimensional  incompressib 

laminar  boundary  layer  flow  over  a  wedge  (i.e.,  Falkner-Skan  flow)  which 

is  subjected  to  a  weak  free  stream  disturbance.  The  disturbance  is  spanwise- 

periodic,  extends  indefinitely  far  downstream,  and  has  a  spanwise  wavelength 

that  is  constant  but  arbitrary.  The  governing  equations  are  formulated 

in  the  orthogonal  coordinate  system  shown  in  Figure  1  with  the  origin  at 

the  leading  edge.  We  introduce  the  spanwise- s inuso idal  disturbance  as 

a  small  perturbation  on  a  basic  2-D  host  flow  U  ,  V  as  follows: 

o  o 


u  (  x ,  y  ,  z  ) 


v  (  x  ,  y ,  z  ) 

w  ( x ,  y ,  z ) 
p(x,z) 


u0U,y) 
vo ( x , y ) 


+  c  u^(x,y)  cos(-^jp) 
+  e  vi  (  x ,  y )  cos(^p) 
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e  w±(x,y)  sin(~ —)  +  ... 
e  p^(x)  cos(—^-)  +  ... 
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where  £<£<1.  Substituting  equation  (3)  into  the  3-D  governing  Navier-St:- 

flow  equations,  neglecting  terms  of  second  order  and  higher  (i.e.,  assuming 
a  weak  disturbance  field),  ar.d  implementing  the  aforementioned  assumptions, 
the  first  order  perturbation  equations  for  the  ,  V^,  and  W^  are: 


J'J:  ovi  2t 

—  +  —  ^~wl  =  0 


-2- 


3u0 

1  ox 

3u0  1 
+  vloy“  +  p 

dPi  Jd2ul 

3x~"  Ldy2 

3w^ 

2  ir  Pi 

P2vl  /2«f  1 

0  dy 

-  7T  T  =  v 

Iay2  "kH 

subject  to  the  boundary  conditions 

U1  =  V1  =  W1  =  o 

U]_  =  0 ,  wi  =  we  @ 


y  =  0 

y  -♦  03 
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Similarly,  the  perturbation  shear  must  vanish  at  the  boundary  layer  edge: 
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Notice  that  the  outer-boundary  conditions  (9)  further  imply  that  (J u^/<)x) 
in  the  inviscid  flow. 


=  0 


We  next  restrict  our  attention  to  the  case  where  Ue  takes  a  power 
law  form  Ue  =  c(x/L)m  and  consequently  the  host  flow  is  of  self-similar 
type  (Falkner-Skan  flow).  Now  this  well-known  solution  is  a  self-similar 
one  in  terms  of  the  stretched  coordinate  "£=  yyCm+l ) 
with  stream  function  and  velocity  components 
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in  which  fQ(n)  is  governed  bv  :: 
split  boundary  value  problem: 


following  ordinary  differential  equation 
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f0(0)  =  f o  <  o )  =  0 
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As  suggested  by  the  foregoing  plus  some  preliminary  study,  we  postulate 
the  solution  for  disturbance  velocity  components  in  the  following  series 
expansion  form: 


(15) 


(  Lo  ) 


Strtiere  is  a  rescaled  non-dimensional  streamwise  distance.  Then 

substituting  equation  (15)  in  (7)  and  setting  the  net  coefficient  of  each 
power  of  f  eq^al  to  zero,  we  obtain  the  following  ordinary 

differential  equations  for  the  cross  flow  functions  G^,  G^. . .  etc: 
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subject  to  the  boundary  conditions 
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Likewise  from  (6)  and  (16) 


the  streamwise  functions 
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with  the  boundary  conditions 
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F  i  ( 0 )  =  f[(0)  =  0  ;  F i ( “ )  -  0  (i-1,2,3,  ...) 


These  disturbance  equations  together  with  the  host  flow  (14)  have  been 
simultaneously  integrated  numerically  by  the  Runge-Kutta  method  [3],  Figures 
2  and  3  show  the  typical  results  for  the  streamwise  and  cross-flow  disturbance 
velocity  functions  F'(^)  and  G respectively ,  with  different  values 
of  m.  Additional  overall  properties  of  these  solutions  are  summarized 
in  Figure  4. 


3.  Discussion  of  Results 


We  observe  from  these  results  that  the  disturbance  velocity  functions 
and  the  corresponding  wall  shears  are  strongly  amplified  by  an  increasingly 
adverse  pressure  gradient.  Furthermore,  we  expect  that  the  first  few  terms 
of  the  series  suffi"  to  give  numerical  solutions  with  desired  accuracy 
because  with  the  in.  casing  order  (i)  of  the  functions  F^,  G^,  Fj_,  and 
Gj_,  their  signs  are  alternating  while  the  magnitudes  of  their  peak  amplitudes 
are  decreasing  for  all  values  of  m-  Thus  we  see  that  an  increasingly  favorable 
pressure  gradient  decreases  the  sensitivity  of  the  boundary  layer  properties 
to  the  3-D  disturbances  whereas  an  adverse  pressure  gradient  dramatically 
amplifies  the  disturbances.  Moreover,  from  Figure  4  and  the  attendant 
skin  friction  relations  (derived  from  the  above  analysis) 
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we  see  thaC  an  increasingly  adverse  pressure  gradient  decreases  the  basic 
term  £q(0)  of  Cf  and,  at  the  same  time,  increases  its  disturbance  component 
Fl(0)  and  so  may  cause  Cf-^0  at  spanwise  stations  z  =^v  [ ( 2n- 1 ) / 2 ]  n=l,2,3,.. 

Ac  these  values  of  z,  Cf  also  vanishes.  Therefore,  we  expect  a  spanwise-per iodi< 
pattern  of  incipient  separation  as  we  intensify  the  adverse  pressure  gradient. 

In  the  full  paper,  the  corresponding  3-D  displacement  thickness  distribution 
will  also  be  examined  in  detail. 

Further  study  of  our  solution  reveals  that  Fannelop's  singularity 
in  which  the  3-D  effect  on  u  grows  downstream  without  bound  is,  in  fact, 
due  to  his  neglect  of  the  perturbation  pressure  and  viscous  z-derivative 
effects,  whereas  the  complete  solution  including  these  terms  taker  a  form 
of  infinite  series  with  Fannelop's  solution  as  only  the  leading  term. 

Thus,  the  perturbation  pressure  and  the  viscous  z-derivatives  correct  and 
bound  the  leading  term  approximation  which  is  otherwise  singular.  Furthermore, 
it  is  seen  that  these  terms  become  more  important  for  problems  with  smaller 
wavelength  disturbances. 
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ABSTRACT 

The  generation  of  streamwise  vorticity  within  the 
laminar  boundary  layer  along  a  spanwise- rippled  plate 
is  theoretically  investigated  for  the  case  where  the  rip¬ 
ple  amplitude  is  small  and  grows  linearly  in  the  stream- 
wise  direction.  Analytical  solutions  are  obtained  for  the 
vortex-induced  disturbance  flow  for  arbitrary  values  of 
the  spanwise  wavelength  to  boundary  layer  thickness 
ratio.  These  are  then  used  to  examine  the  associated 
behavior  of  the  3-D  skin  friction  and  displacement  thick¬ 
ness  distributions. 

NOMENCLATURE 


a(x)  ~  peak  amplitude  of  ripples 

C/.,  Cf,  ~  component  net  skin  friction  coefficients 

fi  =  streamwise  similarity  variables 

y<  =  cross-flow  similarity  variables 

L  =  characteristic  length 

p  a  static  pressure 

R,m ,  R,x  *  Reynolds  numbers  based  on  x  and  A 

u,v,w  at  boundary  layer  velocity  components 

U,  V,  W  si  outer  invisdd  flow  velocity  components 

t,y,t  -  streamwise,  normal  and  cross-flow  coordinates 

SI,  SI  »  component  displacement  thicknessess 

*30  *3-D  displacement  thickness 

<  =  small  expansion  parameter 

Cm,  (»,  C *  —  net  vorticity  components 

q  =  similarity  coordinate 

A  =  wavelength  of  ripples 

I*  *  coefficient  of  viscosity 

v  *  kinetic  viscosity  coefficient 

(  «  nondimensional  streamwise  coordinate 

p  a  density 

r  at  shear  stress 

d  *  harmonic  potential 

4  at  amplitude  of  d 

d  *  stream  function 

subscripts 

e  a  poroperties  at  the  boundary  layer  edge 

w  a  properties  at  the  wall 

oo  a  freestream  quantities 

0  a  undisturbed  flow  properties 

1  a  perturbed  flow  properties 

inv  a  wavy  wall  invisdd  disturbance  solution 
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1.  INTRODUCTION 

The  significant  effect  of  streamwise- vortex  arrays  on 
a  host  boundary  layer  is  well  known.  In  particular,  the 
resulting  enhanced  mixing  rate  of  the  combined  flow 
(“vortex-generator”  effect)  is  of  great  practical  inter¬ 
est  in  conjunction  with  both  external  aerodynamic  flow 
fields  and  those  within  turbomachinery  devices. 

Recently,  Werle  and  his  co-workers  1  have  drawn  at¬ 
tention  to  the  favorable  effects  occurring  from  the  use 
of  a  spanwise- periodic  surface  pattern  on  an  airfoil  to 
generate  such  vortex  arrays  and  have  reported  on  ex¬ 
perimental  observations  of  the  resulting  enhanced  mix¬ 
ing  and  separation  effects  in  the  wake  downstream  of 
the  airfoil.  In  the  present  paper,  we  give  a  theoret¬ 
ical  analysis  of  such  an  array  for  the  simplified  case 
of  attached  incompressible  laminar  boundary  layer  flow 
along  a  slightly  spanwise  rippled  flat  plate  (see  Fig. 
1).  Our  main  focus  is  to  illuminate  the  physics  of  the 
streamwise  vortex  generation  and  its  three-dimensional 
influence  on  the  flow  within  the  highly  viscous  region 
near  the  surface. 

2.  FORMULATION  OF  THE  ANALYSIS 

Our  theoretical  treatment  rests  on  the  following  as¬ 
sumptions: 

s  Small-amplitude  surface  ripples  that  grow  linearly 
with  streamwise  distance:  =  easin  with 

ass  and  <  a  given  non-dimensional  small  param¬ 
eter. 

t  Arbitrary  given  spanwise  wavelength  A 

•  High  Reynolds  number  laminar  flow 

•  Steady  incompressible  flow 

•  Resulting  disturbance  vortices  extend  infinitely  far 
downstream. 

The  analytical  approach  is  a  perturbation  method, 
analogous  to  that  used  by  the  Junior  author  in  a  similar 
study  of  freestream  vortex-array  effects  *,  wherein  the 
flow  properties  are  expressed  as  small  spanwise- periodic 
disturbances  upon  the  basic  Blast  us  boundary-layer: 

=  ->-<Pi(*)sin(^)  +  •••  (1) 

»(*.¥,*)  =  uo(*,y)  +  eui(*,y)sin  +  •••  (2) 

■  t>o(*,y)  ev,(*,y)sin  +•••  (3) 


»(*.».*)  =  «»t(*,y)coe  ("J“) - (liiwMM*)  (4) 

where  ««  and  ««  belong  to  the  Blatius  Row.  Subatitu- 
tion  into  the  basic  governing  Navier-Stokes  equations 
under  the  high  Reynolds  number  conditions  of  negligi¬ 
ble  dp/dy  and  n  (d*  /  0  **)  effects  within  the  boundary 
layer,  retaining  only  the  leading  order  «  approximation 
tor  the  waryn ess- induced  disturbances,  and  then  sub¬ 
tracting  out  the  basic  undisturbed  flow  relations  yield 
appropriate  sets  of  differential  equations  governing  the 
disturbance  distribution  functions  pj,  ultoi,  and  u>i  (see 
2.2  below).  These  are  to  be  solved  subject  to  the  outer 
boundary  conditions  that  u  and  w  at  the  boundary  layer 
edge  approach  values  given  by  their  near-wail  inriscid 
solution  counterparts  to  the  same  rippled  wall  problem, 
while  along  the  impermeable  surface  they  are  subject  to 
the  no  slip  condition,  giving  to  order  c  that 


*i(*,0)  *  0 

(3) 

ut(x,0)  s  -x^(x,0) 

(8) 

t»t(x,0)  2s  0 

(7) 

2.1  Inviscid  Disturbance  Solution 

Irrotational  inviscid  (potential)  flow  past  the  rippled- 
type  of  plate  shown  in  Fi«-  1  exhibits  in  itself  a  non¬ 
trivial  disturbance  solution.  This  is  necessarily  gov¬ 
erned,  in  general,  by  a  harmonic  disturbance  potential 
di(x,y,s)  such  that  a  -  U,  =  ddi/dx  and  w  =  ddi/d» 
vanish  in  the  uniform  mainstream  U9  far  from  the  plate 
while  on  the  surface  satisfying  the  inviscid  imperme¬ 
ability  condition  that  V^x.O)  =  -U9.  In  view  of  Eqs. 
(l)-(4),  di  i*  further  postulated  to  have  the  spanwise- 
periodic  form  *  eff*#(x,y) sin  whereupon  ♦  is 
governed  by  the  Helmholts  equation 


31*  &•*  /2*\» 

dx1  +  dy3  \  A  /  * 


(«) 


subject  to  the  conditions  that  ♦  vanish  for  large  (s’ 
while  jjj(x,0)  *  1. 

Now  some  preliminary  study  reveals  thst  a  solu¬ 
tion  of  Eq.  (8),  having  the  desired  property  that  at 
(and  hence  M/da)  vanish  for  an  assumed  infinite  down¬ 
stream  length  of  ripples.  Is  given  by  the  separation  of 

variables  form* 


From  this,  we  then  obtain  the  inviscid  tidewash  distur¬ 
bance  velocity  at  the  wall  of  the  present  problem  (which 
is  the  outer  boundary  value  for  the  underlying  boundary 
layer  behavior)  as 


*This  solution  excludes  s  wnsfl  ttresawiss  Issdisg  *dg«  rsgioo 
s  <  (i/Jv)  thst  is  also  excluded  by  our  asgltct  of  viscous  x- 
dscivadve  terms  compered  to  thses  o/s-derivstiv*  is  th«  undsrlyiag 
boaxdsry  layer  ragioe. 


=  o)  =  ^■(*•0)  =  -v»  do) 


2.2  Viscous  Boundary  Layer  Disturbance  Field 
The  perturbation  functions  ut,9i  and  wt  of  Eqs.  (1)- 
(4)  within  the  boundary  layer  are  governed  by  the  fol¬ 
lowing  set  of  equations: 


dr 


d*»  /2tr\ 


(ID 


dU|  duj 

u,a7+00'^' 


duo  dtlfl 


id?, 
p  dx 


d*Ui  /2ir\* 

W  '  \TJ 


(12) 


wl  (13) 


dwi  &to\  2ir  pt  fd*Wi  f7v 

“,^_+BodirTT7  =  ,'['y  “It 

These  equations  are  to  be  solved  subject  to  the  wall 
boundary  conditions  given  by  Eqs.  (5)-(?)  with  a  non- 
**ro  $*(x,0)  =  plus  the  outer  inviscid  flow 

matching  conditions  that 


“i(*.  V  —  oo)  =  0  (14) 


wt(r,y  —  oo)  =  1VU..(*,0)  =  -U0  (15) 

corresponding  to  the  vanishing  disturbance  (as  well  as 
basic)  flow  shear  stress  conditions 

^(r,y  -  oo)  =  ^-(r.y  -  oo)  =  0  (16) 

The  attendant  pressure  perturbation  is  determined  by 
observing,  from  conditions  (14)-(16)  plus  Eq.  (13),  that 
the  disturbance  pressure  is  directly  associated  with  the 
cross  flow  viscous  shear  term  —  /i(daw/dr>)  and  equal 
to  the  s  independent  constant  non-dimensional  value 


This  value  is  seen  to  be  very  small  for  high  Reynolds 
number  flows  except  with  small  scale  surface  ripple  wave¬ 
lengths.  In  any  case,  it  docs  not  contribute  to  the  x  • 
momentum  Eq.  (12)  because  dp/dx  =  0. 

Since  the  basic  host  flow  is  a  self-similar  one  (the 
Blasius  solution),  it  proves  convenient  to  reformulate 
the  foregoing  disturbance  flow  problem  in  terms  of  the 
appropriate  similarity  coordinate  i;  =  yiU^/ux)1'1 .  The 
baric  flow  stream  function  vq  (such  that  ««  =  dw<t/'dy, 
«s  =  -Mo /dx)  can  then  be  expressed  in  the  form  Vo  = 
(t'f/’o.x)1/,/o(>7)  where  /0  is  governed  by  Blarius's  well- 
known  ordinary  differential  equation, 

/r+j/o/;  =  o  (is) 

where  ()'  s  d()/dtj  and  /0  satisfies  the  split  boundary 
conditions  /#( 0)  =  /o(0)  =  0,  /#(oo)  —  1  with  /#  —  0. 
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Now,  some  preliminary  study  rereali  that  the  eorre- 
i ponding  perturbation  velocities  in  the  present  problem 
must  behave  for  small  *  as  «i  —  xutF(rf),  w,  —  G(i j) 
where  F  and  G  are  (unctions  of  n  only;  accordingly,  to 
cover  the  entire  mage  of  *  ere  postulate  solutions  in  the 
form  of  the  following  series: 

u,  =  - 1  r  */;(•»)]■  (is) 

w,  =■  -Uo.jtrl9iW  (20) 

%M\ 

where  {  =  *!L  is  a  rescaled  streamwise  distance  in 
terms  of  a  conveniently-chosen  characteristic  length  L 
specified  below,  and  /<(»?)  and  *(ij)  are  governed  by 
suitable  ordinary  differential  equations.  Further,  we  in¬ 
troduce  a  3-D  disturbance  stream  function  such  that 


tr-x  <“> 

(thereby  satisfying  Eq.  (11))  and  postulate  for  it  the 
form 


*  -  k‘/7i(>l)  +  If- 1 l/4(s»)J  (23) 

Then  substituting  expressions  (19)-{23)  into  Equations 
(13)  and  (13)  after  transformation  from  (x,y)  to  «,q), 
equating  to  sera  the  net  coefficients  of  each  power  of  (, 
and  choosing  L  a  (X/2*)1Vo./v  in  order  to  obtain  an 
universal  form  of  the  equations,  we  ultimately  arrive  at 
the  following  set  of  Equations  governing  the  g,  and  /<, 
respectively: 

*+j/o*»0  (24) 

(23) 

9t  +  j/ofs  “  2/£p*  *  j*  (26) 


m  1  .  «  , .  J  f  » 

9i  +  j/o*  ~  ~  l)Ai*  *  fc-i 

/T  +  A/f  +  l/TA-j/i/I-o 

/T  +  j(/s/T  -r  3/0/1)  -  /;/,  +  /;* -  0 

/T  +  \u*n + s/fA)  -  2/,/, 


/r + *  (/o/t + (2s  - 1  )/,v,j  -  (.  - 1  )M 
■  fi-t-fi*- 1 


The  corresponding  boundary  conditions  derive  from  Eqs. 


(5)-(7)  and  (U)-(IS)  and  yield: 

ft(  0)  =  *|(0)  =  0 

(32) 

p[(oo)  —  1,  pwl(oo)  —  0 

(33) 

and 

fU  0)  =  -/o(0) 

(34) 

A(0)  =  /L(0)  =  0 

(35) 

/'(«)  =  0 

(36) 

where  t  =  1,2.3, — 

The  foregoing  equations  constitute  a  system  of  cou¬ 
pled  linear  ordinary  differential  equations  with  split  bound¬ 
ary  conditions  that  can  be  readily  solved  numerically 
by  a  standard  shooting  technique  combined  with  the 
Runge-Kutta  integration  method1.  Before  proceeding 
to  a  presentation  and  discussion  of  the  results,  an  im¬ 
portant  general  feature  of  Eqs.  (24)-(27)  governing  the 

cross  flow  should  be  noted:  the  non- homogeneous  terms 
on  the  right  hand  sides  derive  entirely  from  the  viscous 
p(d* u.w/ffs1)  effect  and  its  associated  small  pressure 
disturbance  (Eq.  17),  and  would  otherwise  be  sero  if 
these  effects  were  neglected  a  priori.  Since  the  boundary 
conditions  on  p,-  for  »  >  2  are  completely  homogeneous, 
this  in  turn  means  that  all  the  p*(q)  would  necessarily 
be  identically  sera  as  well.  Thus,  all  the  terms  in  se¬ 
ries  (20)  that  involve  a  non-sero  power  of  (  physically 
represent  the  entire  cumulative  downstream  influence 
of  this  small  viscous  cross  flow  effect.  Likewise,  all  the 
streamwise  disturbance  functions  /,(q)  for  t  >  3  l.e,  the 
right  hand  rides  of  Eqs.  (30)  and  beyond)  derive  solely 
from  including  the  effects  of  the  /r(fl* u/d*1)  term  and 
would  otherwise  be  sero  if  this  effect  were  neglected. 
In  this  Inter  case,  only  the  first  two  terms  of  the  series 

(19)  would  remain;  these  express  the  leading  approx¬ 
imation  to  the  dual  physical  effects  of  (a)  the  stresun- 
wise  velocity  disturbance  due  to  the  surface  ripple  effect 
via  the  no-slip  condition  at  the  wall  [this  being  —  {1/3 
and  hence  the  most  dominant  effect],  followed  by  (b) 
the  disturbance  within  the  boundary  layer  caused  by 
the  overlying  in  viscid  cross-flow  perturbation  induced 
by  the  ripples,  growing  like  (  and  hence  taking  effect 
further  downstream. 

3.  RESULTS  AND  DISCUSSION 
3.1  Disturbance  Velocities  and  Skin  Friction 
Numerical  results  for  the  first  four  (t  m  1,2,3  and 
4)  disturbance  velocity  functions  &(q)  and  AM  across 
the  boundary  layer  are  presented  in  Figures  2  and  3,  re¬ 
spectively;  the  corresponding  shear  stress  distributions 
p’(q)  and  /*(q)  are  illustrated  in  Figs.  4  and  5.  Be¬ 
yond  i  *  2,  it  is  seen  that  the  maximum  value  of  these 
perturbations  decreases  in  magnitude  and  alternates  in 
sign  with  increasing  «,  implying  that  the  higher  order 
terms  in  (  contribute  to  the  series  solutions  (19)  and 

(20)  in  only  a  small  and  cancelling  way  at  larger  down¬ 
stream  distances.  It  should  be  noted  that  it  becomes 
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increasingly  difficult  to  compute  the  decreasing  values 
of  the  functions  at  larger  «  because  even  a  small  round 
off  error  grows  relatively  large  with  respect  to  the  very 
small  target  value. 

The  special  values  f“( 0)  and  f‘(0),  related  to  the 
stream  wise  and  cross  shear  disturbances  on  the  surface, 
are  tabulated  in  Table  I  up  to  i  *  4.  From  these,  the 
physical  skin  friction  components  can  be  reconstructed 
from  the  present  similarity  series  solution,  as  follows. 
For  the  streamwise  skin  friction  coefficient,  we  have 
Cf,  =  2p(du /9y)m/(pfJX)  which  via  the  similarity  trans¬ 
formation  and  Eq.  (19)  leads  to  the  final  expression 

c,.r“’  =  2/T(o)  -  ^  [(l/v;(o) 

+  £rl/r(0)]«a(Hji)  (37) 

where  /£( 0)  =  .3321  pertains  to  the  undisturbed  flat 
plate  boundary  layer.  The  corresponding  span  wise  skin 
friction  component  Cf,  a  2p(dw/dy)a/(pU^)  likewise 
comes  out  to  be 

C,X'*  =  -  [gr'-VO))]  cos  (^)  (38) 

Several  interesting  conclusions  emerge  from  an  inspec¬ 
tion  of  Eqa.  (37)  and  (38).  First,  since  /,"(0)  is  aero 
while  pi* (0)  is  not  (Table  I),  it  is  seen  that  the  leading- 
term  approximation  for  the  rippled  wall  effect  at  small 
(  does  aai  contribute  to  the  streamwise  shear  stress  Cf, 
but  only  to  the  cross  flow  component;  the  next  term 
(»  *  2)  associated  with  the  invisdd  cross  flow  distur¬ 
bance  produced  by  the  ripples,  however,  does  influence 
Cf,.  Second,  we  can  infer  from  Eq.  (37)  that  the  ripples 
will  hasten  the  onset  of  streamwise  separation  Cf,  —  0 
at  those  spaawise  stations  where  the  perturbation  con¬ 
tribution  on  the  right  hand  side  has  a  nega¬ 

tive  value.  Sines  /J’(O),  p,’(0)  and  fj( 0)  are  all  positive 
(Table  I),  this  will  occur  in  the  leading  approximation 
at  stations  where  (2**/A)  *  3*/2,7e/2, . .  by  Eq.  (38) 
and  Fig.  0,  these  correspond  to  spaawise  locations  that 
are  troughs  in  the  ripples  gad  where  also  tbs  cross  flow 
shear  exactly  vanishes. 


Focusing  on  the  later  streamwise  component  as  the  one 
of  primary  interest  here,  we  thus  get  from  Eqs.  (4)  and 
(20)  that 

(ir)  $rM\  (39) 

In  view  of  the  properties  of  the  functions  g“(f)  shown  in 
Fig.  4,  Eq.  (39)  predicts  as  expected  that  the  streamwise 
vortidty  generated  by  the  surface  ripples  is  a  maximum 
at  the  wall.  Consistent  with  our  original  assumption 
that  the  inviscid  velocity  disturbance  field  caused  by  the 
ripples  is  irrotational,  this  equation  also  correctly  yields 
vanishing  vortidty  outside  the  boundary  layer  since  all 
the  y"(Q)  —  0  aa  q  —  oo.  Finally,  we  note  that  since  the 
/“( 0)  are  positive  for  >  <  2,  the  present  theory  predicts 
in  the  leading  approximation  that  the  largest  stream- 
wise  vortidty  generation  occurs  at  the  lateral  stations 
of  maximum  spaawise  slope  and  cross-flow  shear  stress 
(Fig.  6)  as  one  would  expect  on  physical  grounds. 

3.3  Displacement  Thickness  Distribution 

The  three  dimensional  perturbation  field  due  to  the 
span  wise  ripples  also  alters  the  displacement  thickness 
distribution  along  the  plate  and  hence  the  effective  body 
seen  by  the  invisdd  flow.  Since  this  property  may  be  of 
interest  in  subsequent  studies  of  high  speed  boundary 
layers,  we  conclude  by  examining  it  in  the  present  prob¬ 
lem  to  obtain  some  insight  as  to  the  interactive  "vortex- 
geaerator  effect"  on  the  invisdd  flow. 

Now  the  general  3-D  displacement  thickness  distri¬ 
bution  is  defined  by  the  first  order  partial  dif¬ 

ferential  equation* 

^[^.(«m'O]  +  ^[W.(^-^)]=0  (40) 

where  f*  and  are  the  streamwise  and  cross-flow  dis¬ 
placement  thicknesses  defined  by 

and 


3.3  Streamwise  Vortidty  Generation 

One  of  the  important  features  of  the  present  theory 
is  the  prediction  from  first  prindplcs  of  how  streamwise 
vortidty  (,  is  generated  by  the  spaawise  ripple  effect 
acting  deep  within  the  viscous  boundary  layer.  This 
can  be  determined  from  the  foregoing  analysis  by  devel¬ 
oping  the.  basic  relationship  for  C«  in  terms  of  the  sbove 
similarity-series  solutions. 

Now  it  can  be  shown  from  an  examination  of  the  gen¬ 
eral  Navier-Stokes  equations  that,  consistent  with  the 
high  Reynolds  number  boundary  layer  model  equations 
adopted  in  the  present  study  (specifically,  Eqs.  11-13), 
we  should  also  employ  the  following  boundary  layer- 
approximations  to  the  corresponding  vortidty  compo¬ 
nents:  (,  a  9%/dg,  a  -dm /dt  and  C*  -  -0w/9g. 


Upon  substituting  the  expressions  for  to,  a  and  p,  from 
the  above  analysis,  expanding  the  lower  limits  of  Eqs. 
(41)-(42)  in  a  Taylor  series  ebont  y  =  0,  retaining  only 
the  order  <  3-D  effects,  and  expressing  results  in  terms 
of  oar  similarity  variable  formulation,  we  obtain  the  fol¬ 
lowing  expressions: 


w/*)  »<*•-«  «n 


(t) 


( J is.  i  + 

l  imt 


(»:/*)<’'«.  +  {£  c  ■'*  w} 


(43) 


(44) 
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where  the  various  numerical  values  of  the  constants  a0  s 
Urn,  _„(*>-  fo(.n)\,  =  lim*  -  -  [  H  -  Ji(7)|.  *(<») 

and  /,( 30)  are  given  in  Table  I. 

Now  consistent  with  the  perturbation  approach  of 
the  present  analysis,  the  solution  of  Eq.  (40)  takes  the 
form  i\o  i  fj(*)-«<((x)  iin(2TzA);  after  some  elemen¬ 
tary  calculus,  the  following  results  are  obtained  from 
Eqs.  (43)  and  (44): 

w/*)  (45) 

(Sl'x)R^  =  -<*  -  {g^Bm /,(,) 

+  £  [?*(*)/(*'  *  j)  4-  f -7t(l)]  J  (40) 

It  is  interesting  to  note  that  the  pi  term,  although  con¬ 
tributing  to  S‘,  does  gal  ultimately  appear  in  the  Sf 
expression;  this  is  because  the  lin^_—[j7  - pi(q)j  is  can¬ 
celed  in  Eq.  (40)  by  the  linv-w[t;  -  /o(q)|  term  in  j* 
(see  Table  I).  It  can  be  seen  from  Eq.  (46)  that  the  per¬ 
turbation  of  the  viscous  displacement  effect  due  to  the 
surface  ripples  (in  the  leading  approximation)  it  180* 
out  of  phase  with  the  wall  ripples  and  so  causes  a  local 
increase  in  f  *  in  the  surface  valleys,  while  a  correspond¬ 
ing  thinning  of  the  boundary  layer  occurs  in  the  same 
ripple  pattern  location  as  does  the  minimum  stream  wise 
skin  friction. 


4.  CONCLUDING  REMARKS 

The  present  study  has  shown  that  it  is  possible  to 
construct  a  basic  analysis  of  how  small  amplitude  spanwise- 
periodic  ripples  on  a  surface  generate  streamwise  vor- 
ticity  in  an  overlying  laminar  boundary  layer  flow.  The 
results  should  prove  usefull  as  an  intreptive  guide,  and 
perhaps  also  as  an  upstream  starting  solutions,  in  the 
numerical  treatment  of  the  more  genersd  nonlinear  prob¬ 
lem  associated  with  larger  amplitude  disturbances. 

The  compressible-flow  counterpart  of  the  present  prob¬ 
lem.  including  especially  the  3-D  heat  transfer  distur¬ 
bances  associated  with  these  streamwise  vortices,  would 
be  of  considerable  practical  interest  as  a  follow-on  inves¬ 
tigation. 
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TABLE  I:  Numerical  Results  of  Key  Properties 


»'  m  fl*( 0)  li limn_a.g;(q) 

1  0.0000  0.3321  -0.9954  - 

2  0.1661  0.8756  0.8492  1.090 

3  0.6028  -0.3327  1.558  -0.5782 

4  -0.2343  0.1374  -0.8598  0.08261 
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ABSTRACT 


Tha  treatment  of  turbulence  affacea  oa  super- 
•onic  ■ hock/ turbulent  boundary  layar  interaction 
la  addreaaad  within  tha  context  of  a  trlpla 
daek  approach  valid  for  arbitrary  practical 
Raynolda  numbers  10^  *  Rs^  t  10^®.  Tha  aodaling 
of  tha  addy  vlacoalty  andbaelc  turbulane  boundary 
proli la  affacta  in  aach  dock  la  examined  in  datall 
uaing  Law  of  tha  Wall/Law  of  tha  Waka  concapta  aa 
tha  foundation.  Raaulta  of  paraaatrlc  aeudlaa  ara 
than  glvan,  ahoving  how  aach  of  thaaa  turbulanca 
model  aapacta  lnfluancaa  tha  laportant  interact  ion 
zona  property  of  upstream  influence. 
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apaad  of  aound 

Van  Drieat  wall  turbulanca  damping 
Parana ear  /  . 

akin  friction  coefficient  1 2  UJ  1 

a  ha  pa  factor  4*/e*,  lncoaipraaalbla 
ahapa  factor 

intaractiva  turbulence  affect  on 
inner  deck  thickness 
upatraaa  influence  distance 

Mach  number 

atatlc  praaaura,  atatic  praaaura 
perturbation  (p  -  p_  ) 

•o 

praaaura  Jump  acroaa  undlaturbad 
incident  ahock 

Raynolda  n unbars  based  on  length 
and  boundary  layer  thicknaaa, 

Intaractiva  turbulent  effect  on 
skin  friction  perturbation 
absolute  tas^erature 

basic  intaractiva  vall-turbulence 

peraeeter 

atraanwlaa  and  vertical  Interactive 
disturbance  velocity  consonants 
undisturbed  incoming  boundary  layer 
profile 

streaanrise  and  vertical  coordinates, 
respeclvely 

affective  wall  shift  (dlsplacanant 
height  of  inner  deck)  seen  by 
interactive  invite  id  flow 

J* l 2  l 

specific  heat  ratio 
boundary  layer  thickness 


8  *  boundary  layer  dlsplacsnent  thickness 

8  ,,  inner  deck  sublayer  thicknaaa 
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ordinary  and  kinematic  coefficients 
of  vlacoalty,  reapactivaly 
turbulent  kineaatic  addy  viscosity 
viscosity-temperature  dependence 
exponent  (p'wT*") 
density 

boundary  layer  momentum  thickness 

shear  stress 

interactive  perturbation  of  shear  stress 


Subscripts 


AD 

1 


lnv 

o 

w 


adiabatic  wall  conditions 
undisturbed  lnvlscid  values  ahead  of 
incident  shock 

conditions  at  the  boundary  layer  edge 
lnvlscid  disturbance  soluclon  value 
value  pertaining  to  the  undisturbed 
Incoming  boundary  layer  profile 
conditions  at  the  surface  ("wall") 


1.  INTRODUCTION 

Shock/Boundary  Layer  Interaction  Is  s  signifi¬ 
cant  feature  of  the  flow  fields  around  control 
surfaces  and  within  air-breathing  engine  inlets 
on  suporsonlc/hyper sonic  aerodynamic  vehicles. 
Xxperlmontal  validation  of  CTD  rndes  designed 
to  accurately  predict  the  essential  properties 
of  such  interactions  is  therefore  of  practical 
importance.  In  particular,  it  is  very  desirable 
to  know  tha  sensitivity  of  these  properties 
to  the  basic  parameters  that  govern  the  turbulent 
structure  of  the  Incoming  turbulent  boundary 
layer  (upon  which  the  subsequent  interaction 
depends).  While  the  Influence  of  various  turbulence 
models  oa  different  types  of  aerodynamic  flow 
field  calculations  has  been  extensively  studied, 
especially  under  separated  flow  conditions  (see, 
e.g. ,  Raf.  1),  there  remains  the  need  for  a 
systematic  study  of  the  shock/boundary  layer 
interaction  problem  per  se  within  the  context 
of  the  Lav  of  the  Wall/Wake  concept  that  is 
widely  employed  by  experimentalists  to  characterize 
the  turbulent  boundary  layer.  The  present  paper 
addresses  this  question  for  the  case  of  two- 
dimensional  supersonic  non- separating  shock 
or  compression  corner-generated  interactions 
on  adiabatic  walls. 

Since  it  has  been  clearly  shown  by  the 
late  R.T.  Davis  and  others  that  all  successful 
CFD  treatments  of  such  interactions  must  recognize 
tha  inherent  triple-deck  structure  of  the  inter¬ 
action  zone  (Fig.  1),  our  approach  is  formulated 
in  terms  of  this  structure  using  eddy  viscosity 
concepts.  Because  it  has  proven  applicable 
to  a  very  vide  range  of  Reynolds  numbers  and 
adaptable  to  practical  flow  field  ealcuatlon 
schemes,  we  employ  for  this  purpose  the  non- 
aeymptotlc  version  of  this  triple-deck  theory 
due  to  Inger.2  The  modeling  of  the  eddy  viscosity 
and  mean  velocity  profile  effects  in  each  deck 
is  then  examined  using  the  Lav  of  the  Wall/Wake 
framework  as  the  foundation. 


2.  RATIONALE  Or  THI  T  IIP  LI  DICK  APPROACH 

Since  It  la  tha  foundational  framework  uaad  Co 
addraaa  tha  varloua  curbulenca-aodeling  laauaa, 
a  brlaf  out  Una  of  tha  triple-deck  approach  and  tha 
advantage*  of  lta  naa- aayaptotlc  varalon  will  flrat 
ba  (Ivan.  Wo  conaldar  aaall  dlaturbancaa  of  an 
arbitrary  lncoaing  turbulanc  boundary  layar  dua  to 
a  waak  external  ahock  and  exaaiua  tha  dacallad  par- 
turbaclon  flald  within  tha  layar.  At  high  Raynolda 
nuabara  it  haa  baan  aatabllahad  that  tha  local 

lntaractlon  dlaturbanca  flald  In  tha  aalghborhood 
of  tha  Impinging  ahock  organlxaa  lcaalf  Into 
thraa  baalc  layarad-raglona  or  "deck*"  (Plgura  l)t 

1)  an  outar  raglon  of  potential  lnvlacld  flow 
above  tha  boundary  layar,  which  contains  tha 
Incident  ahock  and  Interactive  wave  ayaceaat 

2)  an  Interaadlace  deck  of  rotatlonal-lnvlacld 
dlaturbanca  flow  occupying  tha  outer  901  or 
aora  of  tha  lncoalng  boundary  layar  thlckneaa t 

3)  an  Inner  aublaycr  adjacanc  Co  the  wall  contain¬ 
ing  both  turbulent  and  laalnar  ahaar  atreaa 
dlaturbancaa,  which  accounta  for  the  Interactive 
akin  friction  perturbatlona  and  hence  any  poaalble 
Incipient  reparation  plua  aoat  of  the  upatraaa 
Influence  of  tha  Interaction.  Tha  "forcing 
fraction"  of  the  problae  hare  la  thua  lnpraaaed 

by  the  outer  deck  upon  tha  boundary  layar i 
tha  aUddla  deck  couplet  thla  to  tha  reeponae 
of  the  Inner  deck  but  In  to  doing  can  ltaelf 
modify  Che  dlaturbanca  flald  to  aoae  extent, 
while  tha  alow  vlacoua  flow  In  tha  thin  Inner 
deck  raacta  vary  atrongly  to  the  praaaure  gradient 
dlaturbancaa  lnpoaed  by  theta  overlying  decka. 

Thla  general  triple  deck  atructura  la  aupported 
by  a  la»e  body  of  experlanntal  and  theoretical 
atudlea.* 

Concarnlng  tha  Importance  of  the  Inner  ahear 
dlaturbanca  deck  and  tha  accuracy  of  deliberately 
ualng  a  non-aayaptocic  Creataent  of  the  detalla 
within  the  boundary  layer,  we  note  that  while 
aayaptotlc  (leg  •*  •  )  theory  pradlcta  an 
exponentlelly-eaall  thlckneaa  and  dlaplaceaent 
affect  contribution  of  tha  Inner  deck,  thla 
la  not  apparently  true  at  ordinary  leynolda 
nuabara,  where  aiany  analytic  and  experimental 
atudlaa  hava  firmly  aatabllahad  that  thla  deck, 
although  Indeed  very  thin,  atlll  contrlbutaa 
algnlf leant ly  to  tha  overlying  Interaction 
and  lta  dlaplacaaent  thlckneaa  growth.*  Thua 
wa  take  tha  point  of  view  hera  that  the  Inner 
deck  la  In  fact  algnlflcant  at  leynolda  nua^era 
of  practical  Interact.  Moreover,  It  contain* 
all  of  the  akin  friction  and  Incipient  aeparatlon 
effecta  In  tha  Interaction,  which  alone  are 
aufflcleat  reaaona  to  exaalne  It  In  detail. 

It  le  further  pointed  out  that  application 
of  aayaptotlc  theory  reaulta  (no  setter  how 
rlgoroua  In  thla  Halt)  to  ordinary  leynolda 
nuabera  la  ltaelf  an  approxlaation  which  aay 
ba  no  aore  accurate  (Indeed  perhapa  laat  ao) 
then  a  phyaleally  well-eonatructed  non-aayuptotlc 
theory.  Direct  extrapolated-aayaptotlc  varaua 
non-aayaptotlc  theory  ceaparlaona  definitely 
ahov  thla  to  be  the  caae  for  laalnar  flova 
(eapeclally  aa  ragarda  the  akin  friction  aapeet) 
and  Che  eltuatlon  can  be  even  woraa  In  turbulent 
flow.  For  exaaple,  the  aayaptotlc  flrat  order 
theory  foraally  excludea  both  the  atreaawl ae 
Interactive  preaaure  gradient  effect  on  the 
ahear  dlaturbanca  deck  and  both  the  normal 


praaaura  gradient  and  ao-callad  "acraaallne  dlver- 
genca"  affacca  on  tha  alddla  dacki  however, 
phyaleal  conalderatlona  plua  experlaantal  obaar- 
vaclona  and  recent  cooperative  numerical  atudlaa  * 
auggaac  chat  thaaa  affacta  are  In  fact  algnlfl¬ 
cant  at  practical  leynolda  nuabara  and  ahould 
not  ba  naglactad.  Of  courca,  aacond  order 
aayaptotlc  corractlona  can  ba  deviaed  to  redreaa 
thla  difficulty  but,  aa  Meyfeh  and  lagab  3 
have  ahown,  run  Che  rlak  of  breaking  down  even 
woraa  whan  extrapolated  to  ordinary  leynolda 
number*.  In  tha  praaanc  work,  w*  avoid  eheae 
problaaa  by  ualng  a  deliberately  nonaayapcoclc 
triple-deck  model  appropriate  to  raallatlc 
Raynolda  nuabara  that  Include*  tha  Inner  deck 
praaaura  gradient  terma  plua  the  alddla  deck 
Hp/ay.  end  atraaallne  divergence  effect*, 
along  with  aoae  aiapllfylng  approximation  that 
render  the  raaultlng  theory  tractlble  froa 
an  engineering  acandpoint. 

3.  TURBULENCE  MODELING  ACROSS  THE  INTERACTION 


Excluding  any  freaatreaa  turbulence,  there 
la  no  explicit  aodellng  needed  in  thla  upper 
raglon  of  potential  lnvlacld  aotloni  the  Influence 
of  the  turbulent  nature  of  the  flow  la  felt 
only  indirectly  through  the  dlaplaceaent  effect 
froa  the  underlying  decka.  The  latter  la  introduced 
by  tha  phyaleal  coupling  condition*  that  both 
v'/U0  and  p*  be  contlnuoua  with  their  alddla 
deck  counterpart*  along  y  9  40- 


Our  analyala  of  thla  layer  reata  on  tha 
key  aiapllfylng  aaauaptlon  that  for  non-aeparatlng 
interaction*  the  turbulent  Reynold*  ahear  atreaa 
change*  are  aaall  and  have  a  negligible  back 
effect  on  the  mean  flow  propertle*  along  the 
Interaction  conei  hence  chit  acre**  can  be 
taken  to  be  "frosan"  along  each  atreealine 
at  It*  appropriate  value  in  the  undlaturbed 
lncoalng  boundary  layer.  Thla  approxlaation, 
llkawlae  adopted  by  a  number  of  earlier  Inveetiga- 
tor*  with  good  reaulta,  la  aupported  not  only 
by  aayaptotlc  analyala  but  eapeclally  by  the 
reaulta  of  lose' a  detailed  experlaantal  atudlea  " 
of  a  non-aeparatlng  ahock  turbulent  boundary 
layer  Interaction  which  thowed  that,  over  the 
ahortranged  interaction  length  atraddllng  the 
ahock,  the  preaaure  gradient  and  Inertial  force* 
outalde  a  thin  layer  near  the  wall  are  at  leaat 
an  order  of  aegnltud*  larger  than  the  correaponding 
change*  In  leynolda  atreaa.  Furtharaore,  there 
la  a  aubatantlal  body  of  related  experimental 
reaulta  on  turbulent  boundary  layer  reapona* 
to  varloua  kind*  of  audden  perturbation*  and 
rapid  preaaure  gradlenta  which  alao  atrongly 
aupport  thla  view?.  Theae  atudlea  unanlaoualy 
confirm  that,  at  leaat  for  non-aeparatlng  flow*, 
algnlflcant  local  Raynolda  ahear  atreaa  dlaturbancec 
are  eaaentlally  confined  to  a  thin  aublayer 
within  the  Law  of  the  Wall  region  (ae*  below) 
where  the  turbulence  rapidly  adjuata  to  the 
local  preaaure  gradient,  while  outalde  thla 
region  where  the  Lav  of  the  Wake  prevail*  the 
turbulent  atreaeoa  reapond  very  alowly  and 
remain  nearly  froten  at  their  Initial  value* 


3.1)  The  Outer  Peck  Flow 


3.2)  Turbulence  Effect*  In  the  Middle  Deck 
Jfr£ien_ Turbulence  Approxlaation 
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far  out  of  tho  local  equilibrium  with  tho  woll  stress. 

Confining  attention,  than,  to  tha  short 
ranga  local  shock  Interaction  tona  where  tha  afore- 
aentloned  "frosan  turbulence"  approximation 
is  applicable,  tho  disturbance  field  caused 
by  a  weak  shock  la  one  of  saall  rotational 
Invlscld  perturbation  of  tha  Incoming  non-unlforn 
turbulent  boundary  layer  profile  M0(y)  governed 
by  tha  aquations 
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as  a  result  of  tha  combined  partlcla-lsantropic 
continuity,  x-mosMntun  and  anargy  eonaarvation 
statements.  It  la  noted  that,  consistent  with 
tha  aaauaed  short  range  character  of  tha  inter* 
action,  tha  streamvlsa  variation  of  tha  undisturbed 
turbulent  boundary  layer  propartlas  that  would 
occur  over  this  range  are  naglactad,  taking 
U0(y),  Po<y>  and  H^y)  to  be  arbitrary  functions 
of  y  only  with  S0,  «0*  and  tv  as  constants. 

Mote  that  Eq.  (3)  la  a  general lxat ion  of  Llght- 
hill'a  well-known  pressure  perturbation  equation 
for  non-unlforai  flows*  which  includes  a  non-linear 
correction  tarsi  for  possible  transonic  effects 
within  tha  boundary  layer  Including  tha  diffracted 
lsiplnglng  shock  above  tha  aonlc  level  of  the 
lncoailng  boundary  layer  profile.  Bqa.  (l)-(3) 
apply  to  a  wide  ranga  of  Incoming  boundary 
layer  profiles  and  provide  an  account  of  lateral 
pressure  gradients  acroas  the  Interactive  boundary 
layer. 


Incoml£g_TJibulen£  iovmdary_Uyr£r_Prof^lle 


The  incoalng  undisturbed  turbulent  boundary 
layer  Is  assumed  to  be  two-dlnanslonal  In  the 
x-dlrectlon  and  to  possess  tha  classical  Lav 
of  the  Wall/Lav  of  the  Hake  structure.  It 
Is  aodaled  by  Vais's*  coaposlta  analytical 
expression  for  tha  resulting  velocity  profile 
coa&lnad  with  an  adiabatic  wall  rafarancs  tempera¬ 
ture  as t hod  correction  for  coapraaslblllty  , 
allowing  arbitrary  non-equlllbrluai  values  of 
Its  shape  factor  Hg.  Thus  If  we  let  *  be  Coles' 
(incoaprasslble)  Hake  function,  n  I  y/A0  and 
denote  for  convenience  I  I  .41  RsjgVfU  ♦  »  ) 
(Tw/T,)l-Hl]  with  •  -  .76  and  Y  -  1.4  for  s 
perfect  gas,  than  the  coaprssalbla  fora  of 
Vais's  coaposlta  profile  nay  be  vrlttem 


2»  ♦  2w  q2. 
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subject  to  the  following  condition  Unking 


»  to  C f  and  Be  *i 

lo  So  ,  _ _ _ 

2»  +  .215  ♦  tn( l+B)  -  .4^/-^(^)  (5) 

Eqa.  (4)  and  (5)  have  the  following  desirable 
propertlasi  (a)  for  n  >.10,  U0/U,  la  dominated 
by  a  Law  of  tha  Wake  behavior  which  correctly 
aatlsfiaa  both  tha  outar  Unit  conditions  UQ/Ue 
and  dll0/dy  ♦  0  as  n  ■  l|  (b)  for  vary  small 
values,  U0  assumes  a  Lav  of  tha  Vail-type  behavior 
consisting  of  a  logarithmic  tarn  that  la  exponen¬ 
tially  damped  out  Into  a  linear  laminar  sublayer 
profile  U/U,  •  ki  as  n  *  0|  (c)  Bq.  (4)  may 
be  diffarantlatad  w.r.t.n  to  yield  an  analytical 
exprasaion  for  dU0/dy  also,  which  proves  advanta¬ 
geous  in  solving  tha  middle  and  lnnar  dack 
Interaction  problama.lt  is  evident  from  these  that 
as  iii,*  1,  the  outer  (wake)  part  of  the  profile 
vanishes  leaving  essentially  a  uniform  (and 
invlseld-llka)  profile  except  for  a  vary  thin 
sublayer  adjacent  to  the  wall. 

Tha  defining  integral  relations  for  <i* 
and  6  i*  yields  tha  following  relationship 
that  links  tha  wake  parameter  to  the  resulting 
compressible  shape  factor  Hi  «  (4i*/8i*)i 


Equations  (4)- (6)  provide  a  vary  general  and 
accurate  modal  of  tha  profile  in  terms  of  three 
Important  physical  quantities i  the  shock  strength 
(Mej.),  tha  displacement  thickness  Reynolds 
number  leg*,  sod  the  Veka  function  s  that  reflects 
the  prior  upstream  history  of  the  incoming 
boundary  layer  including  paatlbl€  aonaquillbrlum 
pressure  gradient  and  surface  mass  transfer 
affects.  Tha  resulting  relationship  of  the 
Incompressible  ahape  factor  Eg  to  the  Vake 
Function  as  a  function  of  Reynolds  number  for 
a  typical  Mi  ■  2.0  flow  is  Illustrated 

in  Fig.  2.  It  Is  seen  from  this  Figure  that 
approaches  a  limiting  value  of  unity  as 
Raj  ♦  «  but  that  this  approach  Is  very  gradual, 
especially  for  vake  function  values  larger 
than  xero  (slightly  favorable  and  adverse  pressure 
gradient  upstream  flow  histories). 

Vlth  these  parameters  prescribed,  the  afore¬ 
mentioned  aquations  may  be  solved  simultaneously 
for  the  attendant  skin  friction  Cf,  tha  value 
of  R  and.  If  desired,  the  Hi  appropriate  to 
thasa  flow  conditions.  Using  tha  adiabatic 
temperature  velocity  relationship 

■  tw,ad  ♦  (T*  *  tw,ad);4^  ^ 

\  !  o, 

the  associated  Mach  number  profile 
Ko(y)  •  U0(YRT0)_li  and  Its  derivative  that 
ars  needed  for  the  middle  deck  Interaction 
solution  nay  then  be  deternln.-d. 

3.3)  Turbulent  Shear  Stress  Disturbances  Alont; 
tha  Inner  Dack 

This  vary  thin  layer  lias  well  within  the 
Law  of  the  Wall  region  of  the  incoming  turbulent 
boundary  layej  profile.  The  original  work 
of  Lighthill1*  treated  It  by  further  neglecting 
the  turbulent  strssses  altogether  and  considaring 
only  tha  laminar  sublayer  effect i  while  this 
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gristly  simplifies  tha  problem  and  ylalda  an 
alagant  analytical  solution,  tha  raaulta  can 
ba  aignifleantly  in  arror  at  high  Reynolds 
nuabort  and  cannot  asp la in  (and  indaad  conflict 
with)  tha  ultiaata  asymptotic  behavior  pertaining 
to  tha  lag  ♦  limit.  Ska  praaant  theory  raaadias 
thia  by  axtaoding  Ughthill's  approach  to  include 
tha  entire  Law  of  the  Hall  region  turbulent 
streea -effects i  tha  resulting  general  shear- 
disturbance  sublayer  theory  provides  a  non- 
asyaptotlc  treatment  which  encoapasses  tha 
complete  range  of  Reynolds  numbers.  It  is 
important  to  note  in  this  connection  that  our 
consideration  of  the  entire  Law  of  the  Hall 
cos^ined  with  tha  use  of  tha  effective  lnviscid 
wall  concept  to  treat  the  inner  deck  displacement 
effect  eliminates  the  need  for  tha  "blending 
layer"  that  is  otherwise  required  to  match 
tha  disturbance  field  in  tha  laminar  sublayer 
region  with  the  middle  lnviscid  decks  except 
for  higher  order  derivative  aspects  of  asymptotic 
matching,  our  inner  solution  effectively  includes 
this  blending  function  since  it  imposes  a  boundary 
condition  of  vanishing  total  (laminar  plus  turbu¬ 
lent)  shear  disturbance  at  the  outer  edge  of 
tha  deck. 

To  facilitate  a  tractlble  theory,  we  Introduce 
the  following  simplifying  assumptions,  (a)  Tha 
incoming  boundary  layer  Law  of  the  Hall  region 
is  characterised  by  a  constant  total  (laminar 
plus  turbulent  eddy)  shear  stress  and  Van  Drieat- 
Cabecl  type  of  damped  eddy  viscosity  modal. 

This  modal  is  known  to  ba  a  good  ona  for  a 
vide  range  of  upstream  non- separating  boundary 
layer  flow  histories,  (b)  For  weak  incident 
shock  strengths,  the  sublayer  disturbance  flow 
is  assumed  to  ba  a  small  perturbation  upon 
tha  incoming  boundary  layer j  la  tha  resulting 
linearised  diaturbanca  equations,  however,  all  the 
physically- Important  effects  of  stresmwlse 
pressure  gradient,  straamvlse  and  vertical 
acceleration,  and  both  laminar  and  turbulent 
disturbances  stresses  are  racaiuadi  (c)  Por 
adiabatic  flown  the  undisturbed  and  perturbation 
flow  Mach  numbers  are  both  quite  small  within 
tha  shear  diaturbanca  sublayer]  consequently, 
the  density  perturbations  in  the  sublayer  dlsburb- 
ance  flow  may  be  neglected  while  the  corresponding 
andast  compressibility  effeet  on  the  Lav  of 
tha  Hall  portion  of  the  undisturbed  profile 
is  quite  adequately  treated  by  tha  Keksrt  reference 
temperature  method  wherein  incompressible  relations 
are  used  based  on  wall  recovery  temperature  proper¬ 
ties  (this  is  equivalent  in  accuracy  to,  but 
easier  than,  the  nan  of  Van  Driest 's  compressible 
Law  of  tha  Hall  profile  7).  (d)  The  turbulent 

fluctuations  and  the  small  interactive  disturbances 
are  assumed  uncorrelated  in  both  the  lower 
and  middle  decks,  (a)  Tha  thinness  of  tha  inner 
deck  allows  the  boundary  layer-type  approximation 
of  neglecting  its  lateral  pressure  gradient. 


where  pw  and  vw  ere  evaluated  at  tha  adiabatic 
wall  recovery  temperature  and  where  It  should 
be  noted  that  tha  kinematic  eddy  viscosity 
perturbation e y  is  being  taken  into  account. 

Tha  corresponding  undisturbed  turbulent  boundary 
layer  Lav  of  the  Hall  profile  U0(y)  is  governed 

by  hi. 


r  (y )  •  const.  •  tw 


d*  (10) 


where  according  to  the  Van  Drlea t-Cebec 1  eddy 
viscosity  model  with  y-t-  *  (y^rWQ/p  Wq)/v Vq 


{.4iv  (1- 


which  yields  for  non- so para ting  flow  disturbances 
that 

..  ,  du 

cT  .  t.41yti-e'y  /  )|2  ~  (12) 
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*T  '  (<JU0/dy'  To 

Here,  A  la  the  so-called  Van  Driest  damping 
"constant]"  wo  use  the  coamonly-accapted  value 
A  ■  26  although  it  is  understood  that  a  larger 
value  nay  improve  tha  experimental  agreement 
la  regions  of  shock- boundary  layer  interaction. 
Substituting  (13)  into  (9)  wa  thus  have  the 
disturbance  momentum  equation 

„  -i,£. 

uo  »x  dy  w  ax 
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from  which  we  have  seen  Chat  inclusion  of  the 
eddy  viscosity  perturbation  has  exactly  doubled 
the  turbulent  shear  stress  disturbance  term. 

He  solve  these  Equations  subject  to  the  wall 
boundary  conditions  Uo(0)  “  u  (x,o)  •  v  (x,  o)  •  0 
plus  an  initial  condition  u  (-«,y)  ■  0  requiring 
that  all  interactive  disturbances  vanish  far 
upstream  of  tha  impinging  shock.  FUthernore, 
at  some  distance  dji,  sufficiently  far  from 
the  vail,  u’  must  pass  over  to  the  lnviscid 
solution  ui'ov  along  the  bottom  of  the  middle 
deck,  as  governed  by 


♦  (.„  )•*  fU  .  o 

dy  v  wo I  dx 


with  Sjl  defined  as  tha  height  where  the  total 
shear  disturbance  (proportional  to  lu'/Jy)  of 
the  inner  solution  vanishes  to  a  desired  accuracy. 

4.  SOLUTION  METHODOLOGY  AND  RESULTS 

The  solution  to  the  foregoing  triple  deck 
problem  is  achieved  for  small  linearised  disturb¬ 
ances  shead  of,  behind  and  below  the  local 
shock  jump,  which  gives  ressonably  accurate 
predictions  for  all  the  properties  of  engineering 
interest.  The  resulting  equations  can  be  solved  by 
a  Fourier  transform  method  to  yield  the  viscous 
interaction  field  physics  for  non- 


••paratlog  flow*  Including  tba  upatreaa  Influence, 
tha  lateral  praaaura  gradient  near  the  a hock 
and  tha  onaat  of  lnclplaat  aaparatlon  (aaa 
la  fa.  2  aad  12  for  tba  datalla  of  thla  aolutloa). 
Nuamrous  detailed  eeaparlsons  with  experiment* 
have  abown  that  It  gives  a  good  account  of  all 
tha  Important  fsatnamo  of  tha  Interaction  over 
a  wide  range  of  hash ‘Reynolds  nuabar  condltlona. 

♦.1)  Fourier  Trane format Ion  Method 

We  only  briefly  outline  hare  tha  atepa  Involve^ 
•Inca  full  datalla  can  be  found  alaawhara.  follow¬ 
ing  Fourier  Trans foraat Ion  w.r.t.x  ,  tha  resulting 
elddla  deck  praaaura  problea  froa  Eq.  3  la  an 
ordinary  differential  equation  In  y  that  can 
be  solved  nuaarlcally  quite  efficiently  for  tha 
Input  turbulent  boundary  layer  profile  Mo(y) 
of  Section  3.2.  In  particular,  for  tha  upatreaa 
Interactive  pressure  rise  we  find  froa  tha 
appropriate  Fourier  Inversion  procasa  using 
tha  calculus  of  residues  that 

■  ip  e*/tu  (16 

where  AP  la  the  overall  shock  pressure  juap 
while  tu  la  tha  characteristic  upatreaa  distance 
given  by 
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in  taraa  of  the  following  profile-dependent 
Integrals  evaluated  by  the  aforeaantlonad  turbulent 
Law  of  tha  Wall/Law  of  tba  Hake  nodal  i 
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The  paraaater  Fweff  here  la  the  effective  lnvtscld 
wall  shift  given  by  the  displacement  thickness 
of  the  underlying  Inner  deck. 

The  corresponding  Fourier  trans foraat Ion 
of  the  Inner  dock  problea  of  Section  3.3,  followed 
by  the  Introduction  of  new  Inner  deck  variables  and 
y-scallng  defined  by  Inger*,  yields  a  set  of 
ordinary  differential  equation  boundary  value 
problems  In  a  "universal"  fora  that  can  be  solved 
and  tabulated  once  and  for  all.  An  example  of  this 
la  Illustrated  In  Fig.  3,  which  shows  the  resulting 
Inner  deck  streaawlse  velocity  profiles  la  tsras 
of  the  eddy  vlscoststy  effect  as  expressed  by  tha 
authors'  Interactive  Turbulence  Parameter2 
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The  typical  laynolda  nuabar  and  waka 

function-dependence  of  this  paraaater  Is  Illustrated 
la  Figure  4,  where  It  la  seen  that  it  grows  to  large 
values  with  Increasing  ,  as  well  as  Increasing 
with  a  . 

we  further  obtain  the  following  result  for 
the  deck's  dicplacesMat  thickness i 


where  tha  eddy  viscosity  effect-function  H(t) 
is  given  In  Figure  3.  The  simultaneous  solution 
of  Eas.  (17)-(19)  for  iu  and  Fveff  Implements 
the  matching  of  tha  Inner  and  middle  decks. 

Tha  resulting  values  of  the  Inner  deck  height 
expressed  as  a  fraction  of  tha  Incoming  undisturbed 
boundary  layer  thickness  are  plotted  versus 
laynolds  nuabar  with  s  as  a  parameter  in  Figure  6 ; 
also  shown  for  comparison  are  tha  corresponding 
sonic  height  ratio  values.  It  is  clearly  seen 
how  rapidly  Pweff^o  decreases  with  Increasing  Re  , 
reaching  exceedingly  saall  values  Indeed,  relative 
to  the  much  more  gradual  decrease  In  y>0nlc/<o- 
It  Is  also  Interesting  to  note  hare,  as  one 
would  expect  on  physical  grounds,  that  while 
the  Inner  deck  thickness  la  hardly  affected 
by  s,  tha  sonic  height  (which  lies  within  the 
wake  region)  ,1s,  significantly  Influenced  and 
increases  with  the  value  of  the  Wake  function. 

Finally,  we  note  the  companion  result  for 
the  upstreaa  skin  friction  that 


w'  •  -1.372  p„'  (x)f  (») 

P  i 

hare 
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and  S(T)  Is  another  Interactive-turbulence 
effect  function,  also  plotted  In  Figure  5. 

Figure  5  Is  a  central  result  of  the  present 
general  turbulent  shear-disturbance  Inner  deck 
treatment;  It  gives  a  unified  account  of  the 
Inner  Interactive  physics  over  the  entire  Reynolds 
nuabar  range  from  quasi- laalnar  behavior  at 
T  <<  1  (lower  Reynolds  nuabers)  to  tha  opposite 
•stream  of  wall  turbulence-doalnatsd  behavior 
atT  »>  1  pertaining  to  asymptotic  theory  at 
very  large  Reynolds  numbers  where  the  inner 
deck  thickness  and  Its  disturbance  field  become 
vanishingly  small. 

4,2)  Predictive  Results  Showlna  the  Role 

of  tha  Turbulence  Mods  ling  Parameters 

A  computer  program  has  been  constructed 
to  carry  out  tha  foregoing  solution  method, 

It  Involves  the  middle-deck  disturbance  pressure 
solution  coupled  to  the  inner  deck  by  means 
of  the  affective  wall  shift  combined  with  an 
upatreaa  Influence  solution  subroutine  (the 
corresponding  local  total  Interactive  displacement 
thickness  growth  and  skin  friction  are  also 
obtained).  This  provides  a  very  general  fundamen¬ 
tal  description  of  the  boundary  layer  In  terms 
of  three  arbitrary  parameters ■  preshock  Mach 
nuabar,  boundary  layer  displacement  thickness 
Reynolds  nuabar,  and  either  the  waks  function 
ror  the  Incompressible  shape  factor  81} • 

Rased  on  the  aforementioned  prograa,  an 
extensive  parametric  study  has  been  carried 
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ouc  to  show  eh«  sensitivity  of  predicted  lnterec- 
C loo  ion*  propareloa  to  eh«  various  key  turbulent 
flow  modeling  parameters.  for  example,  Figure  7 
shows  for  a  typical  ■  2.0  Interaction  the 
affact  of  the  Waka  Function  on  the  interactive 
up* trees  influence  distance  tu  (in  ratio  to 
the  undisturbed  boundary  layer  thickness  4„) 
as  a  function  of  Reynolds  number.  Clearly, 
the  effect  la  as  important  one  over  a  wide  range 
of  laynolda  numbers,  indicating  a  significant 
Increase  in  tu/50  with  r,  and  suggests  that 
accounting  for  the  upstream  boundary  layer  history 
can  be  important  in  the  experimental  validation 
of  CFO  predictions  of  such  interactions.  In 
this  connection,  it  should  be  noted  that  this 
wake  function  aspect  is  totally  lost  in  the 
leading  approxlMtlon  of  the  asymptotic  triple 
deck  approach  (which  is  based  on  Che  limiting 
value  Hi}  *  1.0  pertaining  to  Che  Infinite  laynolda 
number  limit,  wherein  the  wake  component  completely 
vanishes ) . 

Another  Interesting  aspect  of  the  turbulence 
modeling  is  Che  eddy  viscosity  perturbation 
effect  la  the  lnnter  deck)  this  is  illustrated 
la  Figure  8,  where  we  show  how  the  predicted 
upstream  lnfluance  distance  is  altered  by  includlte 
(or  neglecting)  this  affect.  At  moderacely-hlgh 
laynolda  mnbers  (la,  i  10*),  Che  effect  la 
seen  to  be  quite  large,  such  that  neglect  of 
the  interactive  disturbance  to  ey  can  consequently 
underpradlct  iu  by  hundreda  of  percent .  On 
Che  other  hand,  at  very  large  la  where  the  inter¬ 
active  flow  is  essentially  lnviacid-domlnatad 
and  Influenced  only  by  the  outer  waka  region 
of  the  incoming  boundary  layer,  the  eddy  viscosity 
perturbations  have  only  a  small  effect.  Figure  8 
also  serves  to  rees^healse  the  fact  that  the 
present  theory  applies  to  a  very  wide  range 
of  practical  laynolda  numbers. 

In  Figure  9,  we  conclude  by  illustrating 
the  excellent  agreement  of  the  preaent  predictions 
of  tu/<o  with  experiment  When  one  properly 
accounts  for  the  important  effect  of  the  Incoming 
boundary  layer  shape  factor  (or  waka  function). 
Clearly,  the  experimental  validation  of  any 
theoretical  prediction  will  require  a  rather 
careful  determination  of  the  wake  component 
aspects  of  the  incoming  turbulent  boundary  layer. 
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Fig.  4.  Variation  of  tha  Interactive 
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Reynold*  Number  for  Various  Wake 
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Fig.  3.  Streamwlae  Disturbance  Velocity 
Profiles  Across  the  Inner  Deck 
for  Various  Values  of  tha 
Interactive  Turbulence  Parameter 


Fig.  5.  Turbulent  Interaction  Parasntar 
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Fig.  6.  Non-Dimensional land  Im*r  Deck 
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SUMMARY 

RESEARCH  NEEDS  IN  3-D  SHOCK-WAVE/ 

TURBULENT  BOUNDARY  LAYERS  INTERACTIONS 

After  considerable  discussion,  thirteen  separate  topics  were  identified  as 

requiring  immediate  research  attention.  They  are  (in  no  particular  order); as 

1.  Analytical  treatments  of  near-wall  boundary  layer  behavior  to  develop 
efficient  wall -functions  for  computation. 

2.  Experimental  (and  computational)  studies  to  determine  the  effect  of 
modifying  the  upstream  boundary  layer  using  pressure  gradients,  blowing 
and  roughness. 

3.  Measurements  of  skin  friction  and  heat  transfer  distributions. 

4.  Studies  to  identify  where  turbulence  models  are  important  for  the 
accuracy  of  the  computations  (in  large  regions  of  the  interaction, 
turbulence  may  not  be  playing  a  significant  role) . 

5.  Extension  of  calculations  and  experiments  to  higher  Mach  numbers 

(  >  6  ). 

6.  Extension  of  calculations  and  experiments  to  more  complex  interactions 
(for  example,  shock/shock  interactions,  interactions  with  floor  and 
sidewall  boundary  layers,  more  complicated  shock-generator  geometry). 

7.  More  detailed  flowfield  investigations  for  carefully  selected 
interactions.  These  measurements  should  have  enough  detail  to  give 
the  flowfield  structure,  unsteady  characteristics,  turbulence  behavior, 
skin  friction  and  heat  transfer  distributions. 

8.  Studies  of  the  downstream  boundary  layer  relaxation  (very  important  for 
inlet  design) . 

9.  Alteration  and  control  applied  within  the  interaction  zone  using,  for 
example,  bleeding,  blowing,  flow  guides. 

10.  Exploratory  calculation  of  non-adiabatic  interactions. 

11.  Full  characterization  of  incoming  flow  conditions,  including  freestream 
turbulence,  spanwise  inhomogeneities  (3-D  effects  in  nominally  2-D 
flows) ,  changes  in  boundary  layer  turbulence  due  to  Mach  number  effects 
(especially  M  >  6) . 


